Decay of the classical Loschmidt echo in integrable systems 
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We study both analytically and numerically the decay of fidelity of classical motion for integrable 
systems. We find that the decay can exhibit two qualitatively different behaviors, namely an algebraic 
decay, that is due to the perturbation of the shape of the tori, or a ballistic decay, that is associated 
with perturbing the frequencies of the tori. The type of decay depends on initial conditions and on 
the shape of the perturbation but, for small enough perturbations, not on its size. We demonstrate 
numerically this general behavior for the cases of the twist map, the rectangular billiard, and the 
kicked rotor in the almost integrable regime. 
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I. INTRODUCTION 

The study of the relation between classical and quan- 
tum dynamical chaos has greatly improved our under- 
standing of the behavior of quantum systems. An issue 
which may be of interest in several situations is the sta- 
bility of motion. In this respect, even though Liouville 
equation, which governs the evolution of classical dis- 
tribution functions, is linear, the exponential sensitivity 
of classical trajectories with respect to perturbing the 
initial conditions leads to a strong dynamical instabil- 
ity, and characterizes classical chaos. Like the classical 
Liouville equation, the Schrodinger equation is also lin- 
ear. However, the quantum evolution of states is stable 
and this qualitative difference is clearly apparent in the 
Loschmidt echo numerical experiments of Ref. [1]. The 
problem of the stability of quantum motion under per- 
turbations in the Hamiltonian has recently gained a re- 
newed interest [2-16], also in connection with quantum 
computation [17-19]. The quantity of central interest 
in these investigations is the fidelity f q (t) (also called 
Loschmidt echo), which measures the accuracy to which 
a quantum state can be recovered by inverting, at time t, 
the dynamics with a perturbed Hamiltonian. The main 
interest has been focused on classically chaotic systems 
for which, besides numerical experiments, some theoret- 
ical tools are available, like random matrix theory and 
semiclassical methods. However, in the general case, the 
phase space structure is mixed with chaotic components 
and islands of stability. If the motion starts inside an 
integrable island, then it very much resembles the mo- 
tion in integrable systems. Contrary to chaotic systems, 
which are dynamically unstable but structurally stable, 
integrable systems are dynamically stable but very sensi- 
tive to external perturbations. Therefore the analysis of 
the fidelity requires particular care and one may expect it 
to be dependent on initial conditions and on the type of 
perturbation. Indeed, the decay of fidelity in integrable 
systems has been discussed in recent papers [11,20,21], 



and very different behaviors have been found. Jacquod 
et al. have shown the existence of a regime in which the 
quantum fidelity for classically integrable systems decays 
as a power law, with an anomalous exponent of purely 
quantum origin [21]. Prosen and Znidaric have instead 
discussed a regime in which quantum fidelity exhibits a 
much faster Gaussian decay [11]. Both regimes have been 
also discussed by Eckhardt in his analysis of the decay of 
classical fidelity [20] , in which the problem of the evolu- 
tion of classical phase space densities has been addressed 
for linearized flows. 

In the present paper, we discuss the behavior of fidelity 
for integrable classical systems. Besides being of interest 
on its own, our classical study will allow us to understand 
the main mechanisms for the fidelity decay and therefore 
will constitute a valuable reference point for the quantum 
analysis. Here we show the existence of a critical bor- 
der depending on the shape of the perturbation, which 
separates two different types of fidelity decay: a power 
law decay oc l/t n , where n is the dimension of the sys- 
tem, and a much faster decay of ballistic type. We stress 
that the type of decay depends on initial conditions and 
on the shape of the perturbation but, for small enough 
perturbations, not on its strength. We derive an analyt- 
ical expression for the critical border and our theoretical 
results are confirmed by a numerical analysis on three 
different models: the twist map, the rectangular billiard 
and the kicked rotor. 

The outline of the paper is as follows: In Section II we 
develop a general theory for the decay of classical fidelity 
in integrable systems. Section III demonstrates numer- 
ically the validity of our theory in two different typical 
examples of integrable systems, the twist map and the 
rectangular billiard, and in an almost integrable system, 
the kicked rotor. Our conclusions are drawn in Section 
IV. Finally, in the Appendix A, we discuss the long-time 
relaxation to equilibrium. 
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II. THEORY 

The quantum fidelity is denned as the overlap at time 
t of the states \ip(t)) and \ip e (t)), obtained by the evolu- 
tion of the same initial state \ip(0)) with the unperturbed 
Hamiltonian H and the perturbed Hamiltonian H +eV 7 
respectively. The fidelity is then given by 



f g {t) = \W)\Mt))\' 



(i) 



This expression can be equivalently rewritten in terms of 
the Wigner functions as 

/,(*) = (2vh) n J d n q d n p W e (q,p;t) W(a,p;t), (2) 

where n is the number of degrees of freedom. Since the 
Wigner functions can be considered as the quantum ana- 
logues of the classical phase space densities, we define the 
classical fidelity as 



/(*) = J d n qd n p p e (q,p;t) p(q,p;t), 



(3) 



where p, p e are the square normalized classical phase 
space densities (/ d n q d n p p 2 = J d n q d n p p\ — 1). We 
note that f(t) is the classical limit of f q (t). As the den- 
sity evolution is unitary in both classical and quantum 
mechanics, instead of evolving two densities forward in 
time and calculating their overlap, we may first evolve the 
initial density po forward in time with the unperturbed 
Hamiltonian Hq, and then evolve this density backward 
in time with the perturbed Hamiltonian Hq + eV. We 
denote the density obtained in such a way as p 2 t- The 
fidelity is then given by the overlap of the density p 2t 
with the initial density p n : 



/(*) 



= J cTq d n p 



P2t(q,p) po(q,p)- 



(4) 



Such an approach is more convenient for our discussion 
of the classical fidelity of integrable systems. 

For the following discussion we assume that the per- 
turbation of the integrable system is of the KAM type, 
namely that, for small enough perturbations eV , most of 
the tori of the system are only slightly deformed but not 
destroyed. Therefore for most of the tori the transfor- 
mation from old action-angle variables I, to new ones 
I', 0' is possible, in such a way that the new actions are 
constants of motion of the perturbed system. To the first 
order in the perturbation strength e the transformation 
can be written as 



0' = + ef(I,0), 
I' = I + eg(I,0). 



(5) 
(6) 



After the forward unperturbed evolution up to time t 
we have 



& t = ©o + n(i)t. 



(7) 



Then we perform a backward evolution of the perturbed 
system from time t to time 2t, getting 



&' 2t = &' t - n'(i')t. 



(8) 



The frequency vectors il and characterize the linear 
(in time) evolution of the angle variables in integrable 
systems. Expressed in the original action-angle variables, 
the overall evolution can be written as 



24 = O + («(I) - «'(!')) t + 0(e), 



(9) 



where the error term 0(e) is due to the change from one 
set of variables to the other at time t and the reverse 
change at time 2t. Since the perturbation is small, we 
may write the change in frequency as 

n'(i') = n(i) + Afi(i) + ^(i' - 1) + 0(e 2 ), (io) 

where Afi = fi'(I) — O(I) denotes the change of fre- 
quencies on the unperturbed torus I and I' — I gives the 
change of the action variables caused by the perturbation 
(written in Eq. (6) to the first order in e). 

If we consider the angle variables to be uniformly dis- 
tributed at the time t at which the motion is inverted, 
we may introduce the distribution 



Wi 



t) 



= tttV / dn& 
(2tt)« J 



AI -(!'(!, ©)-!)) 



(11) 



which gives the probability density for the transition from 
the torus characterized by the action variables I in the 
unperturbed coordinates to the torus with action vari- 
ables I' in the perturbed coordinates. The e scaling has 
been chosen in order for the function W\ itself not to 
depend on e in the linear approximation. 

In the generic case, the motion on a torus is ergodic. 
It is, however, not random and therefore, in order for 
the equation (11) to well describe transitions between 
the tori, one needs to consider an ensemble of tori in the 
vicinity of the chosen actions I, as only for an ensemble 
of tori with different frequencies we can expect the angle 
variables to be uniformly distributed after a sufficiently 
long time t. Indeed, the spread of frequencies given by 
Sfl = ^fSI translates into the spread of angle variables 
<50 = Sfl t. The time for this spread in the angle variable 
to become comparable to 2ir is 



te 



2n 



\W) Vl 



(12) 



where z/j is the characteristic width of the initial phase 
space density distribution p along the action direction 
(to simplify writing, we have given Eq. (12) for the one- 
dimcnsional case). Thus, our theory based on Eq. (11) is 
valid for times t > tg. 
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Since the final angle variables, after the forward and 
backward evolutions, depend on the actions of the per- 
turbed system, we compute the distribution of the angle 
variables from the distribution of the perturbed actions 



as 



P I (0 2t -0 o ;i) = Wi 



I' -I 



0((I'-I)/e) 



d(& 2t - © ) 



(13) 



Using the expressions (9) and (10), we obtain 

-l 



Pi(02t-©o;i) 



1 



(ct)« 



an 



dl 



X 



Wi 



In 



& -& 2t + O(e) A(l s 
et e , 



(14) 



This expression is the kernel for the combined forward 
and backward evolution of the phase space densities. 

We assume that the width v\ along the action direction 
of the initial density po is much larger than the change 
of the action variable induced by the perturbation, that 
is 



vi > e max Iff (i", 0)1 



(15) 



(to simplify writing we have given the condition (15) for 
the one-dimensional case). Therefore the effects of the 
forward and backward evolutions are felt mainly in the 
change of the angles variable. This means that the evo- 
lution from the initial phase space density po to p 2t is 
given, up to corrections of order e, by 



P2t(l,@) 



J d U & Pi(©' 



&;t) p (I,&'). (16) 



The fidelity f(t) can then be computed by inserting p 2t 
into Eq.(4). 

The kernel Pi(0' — 0) is stretched linearly in time, 
while at the same time it moves ballistically (linearly with 
time) with velocity Att. Under the assumption that the 
perturbation of the shape of the tori is not divergent (as 
it is the case for most of the tori in a KAM regime), 
the distribution Wj(I/e) has a bounded support which is 
determined by the change of the shape of the tori due 
to the perturbation. We can see from Eq. (14) that at 
long times the argument of the function W\ is given by 
— [dfl / dl] -1 (Att / e) . Therefore the long time behavior of 
Pi depends on whether the value of — [dfl/dI]~ 1 (Afl/e) 
falls within the support of W\ or not. In the first case, 
Wi = WiHan/afl-^An/e)) is different from zero and 
therefore the kernel Pi drops cx 1/t™. In the latter case, 
W\ = 0, and therefore Pi drops ballistically. The tran- 
sition between these two regimes is determined by the 
equality 



e 



an 



1 An 



(17) 



where AI s /e are the coordinates of the border of the sup- 
port of Wi. 

We can therefore draw the following conclusions: If 
the perturbation of a classical integrable system is such 
that the primary effect is the change of the shape of the 
tori, then the expected decay of fidelity is cx 1/t™. On 
the contrary, if the change of the frequencies of the tori 
is the dominant effect, then we expect a ballistic decay 
of fidelity, that is the center of mass motion of the phase 
space densities after the forward and backward evolutions 
is responsible for a drastic drop of fidelity. Such a decay 
takes place as soon as the centers of mass of the densities 
pit and po are separated in the angle variables by more 
than their characteristic width As it can be seen 
from Eq. (4), the exact form of the fidelity drop in the 
ballistic regime depends on the tails of the initial distri- 
bution po: for instance, a Gaussian tail gives a Gaussian 
decay of fidelity, whereas a sharp border induces a sharp 
drop to zero of fidelity. Finally, it is important to stress 
that the type of decay, power law or ballistic, depends on 
initial conditions and on the shape of the perturbation. 
However, it does not depend on the strength of the per- 
turbation, provided that it is sufficiently small. Indeed, 
Eq. (17) shows that AI s /e is e-independent (to the first 
order in e), since Aft cx e. 



III. NUMERICAL DEMONSTRATION 

As a first example we consider the perturbed twist 
map, defined by 

I t+ i = It + ecos(a)sin(e t ), 
©t+i =e t + / t+ i+esin(a)sin(/ t+1 ), (18) 

where the angle a determines the mixture between purely 
perturbing the shape of the tori (a — 0) or purely chang- 
ing their frequencies (a = 7r/2). This paramctrization 
allows us to change the type of the perturbation without 
changing its overall magnitude. The change of frequency 
associated with the perturbation is given by 



Afi = esin(a) sin(J). 



(19) 



The conserved action variable of the e-perturbed system 
is, to the first order (in e) approximation, given by 



I' = I + e cos(a) 



(20) 



Indeed, inserting this expression into the mapping (18), 
one can easily verify that It+i = It- The transition prob- 
ability function Wj for this system can thus be obtained 
by means of Eq. (11): 



2tt/ 



dO S 



Wj (Al/e) = 
AI cos(a) 



2sin(7/2) 



cos — 



(21) 
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which gives 



(22) 



cos(a) 
2 sin (7/2) 



(A7/ e ) 5 



The range of the support of the distribution Wj is be- 
tween ±cos(a)/(2sin(7/2)). The critical value a c for 
which AI s /e = — (dCl/dl) _1 Af2/e is therefore deter- 
mined by 



tan(a c ) 



1 



2sin(7) sin(7/2) 



(23) 



In Fig. 1 we show the numerically computed behav- 
ior of fidelity for this system as a function of time for 
various values of the parameter a. We take as initial 
phase space density a rectangle centered around the point 
(9 = 7r, I = 1) with sides of length v® = 2 x 10~ 3 , vi = 
2 x 10~ 2 . The perturbation strength is e = 1CT 6 . To 
compute fidelity, we follow the evolution of N = 10 4 tra- 
jectories, which at t = are uniformly distributed inside 
the above rectangle. The fidelity is then given by the 
percentage of trajectories that return back to this region 
after the forward and backward evolutions. In all cases 
we observe an initial plateau during which the fidelity 
does not decay appreciably. This plateau persists until 
the time t p at which the width of the kernel (14) or the 
shift of its center become comparable to the width v@ 
of the phase space density along the angle variable. In 
cither case this time is 



t p oc 



(24) 



According to Eq. (17), we expect the behavior to change 
from algebraic decay to ballistic one at the value of the 
parameter a = a c . For the chosen initial conditions, 
Eq. (23) gives a c w 0.892. Indeed, the change from an 
algebraic fidelity decay f(t) oc 1/t when a < a c to a 
sharp drop of fidelity when a > a c is clearly seen in 
Fig. 1. 

An interesting feature is that, approaching the criti- 
cal value a c , we observe that the fidelity decay, power 
law or ballistic, sets in after longer and longer times. 
This fact has a clear explanation: The value W\ — 
Wj(— [dfl/dI]~ 1 (Afl/e)), which determines the long 
time behavior of the evolution kernel (14), diverges close 
to the critical value a — a c (see Eq. (22)). When the 
fidelity decay is power law, we have Pi = c/t n , with the 
constant c oc W\. Since c becomes larger and larger close 
to the critical point, the fidelity decay must be postponed 
to longer and longer times. On the other hand, when 
the long time fidelity decay is ballistic, we can see from 
Eq. (14) that the argument of the function W\ goes out- 
side the support of W\ after a time that becomes longer 
close to the critical point. Only after this time the fidelity 
drops off. Of course the decay cannot be postponed in- 
definitely since the exact condition (23) can be satisfied 



only for a single torus, while we always deal with a fam- 
ily of tori upon which the initial phase space density po 
rests. 




FIG. 1. Fidelity decay for the twist map at various val- 
ues of the parameter a — (full line), 0.8 (dashed), 0.892 
(dot-dashed), 1.0 (dot-dot-dashed) and tt/2 (dotted). The 
oc 1/t decay is shown as a thin dotted line 

We also checked numerically that, provided that the 
perturbation is much smaller than the characteristic 
widths vi of the initial density (that is, the requirement 
(15) is fulfilled), the type of behavior does not alter with 
changing the actual size of the perturbation e, as it is 
expected from our theory. We simply rescale the time 
t p oc 1/e after which the fidelity decay starts, in agree- 
ment with Eq. (24). 

To illustrate the fidelity decay in integrable systems 
with more than one degree of freedom, we consider the 
following system: 

h(i u i 2 , e 1; e 2 ) = H (h,i 2 ) + ev(i u i 2 , ex, e 2 ), (25) 



where the unperturbed Hamiltonian 



Ho 



a 2 



*0 



(26) 



describes the motion of a particle bouncing elastically in- 
side a rectangular billiard and the perturbation is given 

by 

V = cos(f3) cos(Oi) cos(9 2 ) + sin(/3)7i/ 2 . (27) 

Again, depending on the value of the parameter /3, the 
perturbation mainly affects either the shape of the tori 
or their frequencies. We use the first order perturbation 
theory of Hamiltonian systems (see, e.g., Ref. [22]) to de- 
termine the effects of the perturbation. What we need 
to find is a set of action-angle coordinates such that, to 
the first order in the perturbation strength e, the Hamil- 
tonian (25) in these new coordinates can be written as a 
function of only the new actions, namely 



H(h, I 2 , 0i, 6 2 ) = 77' (7{, V 2 ) + 0(e 2 ). 



(28) 
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Introducing the generating function 

G(i[,i^,o 1 ,e 2 ) = i[Q 1 + i 2 Q 2 

ecos/3 



sin(6i + 6 2 _) sin(6i - 9 2 ) 



we get 



r _ dG _ 

Il -w 1 - h 



ct\I[ + a 2 I 2 "i^i - a zl 2 
e cos(P) 



(29) 



1 



h = 



8G 
W 2 



a\l\ — a 2 I 2 
ecos(P) 



1 

aih + a 2 I 2 



cos(6i - 6 2 ) 



a\I\ + a 2 I 2 



a\I\ — a 2 I 2 



cos(6i - 6 2 ) 



cos(9i + 6 2 ) 
(30) 

cos(9i + e 2 ) 
(31) 



Substituting the above expressions into the Hamiltonian 
(25), we get 



HVi, 4) = ^I[ 2 + ftf + esin(/3)/i/ 2 . 

The new frequencies are then given by 
B H' 

0' 1 = — =a 1 7{ + esin(/?)7 2 , 



B H 1 

n' 2 = — =a 2 I 2 + esm(P)r i . 



(32) 

(33) 
(34) 



Thus the frequencies changes read as follows: 

Afii = esin(/3)7 2 , (35) 



Afi 2 = esint/?)/!. 



(36) 



As in the previous example, the above expressions allow 
us to find the transition probability function 

Wi(A/i/e,A/ 2 /e) = 
2 1 



71 M ".SOT V / 

y ^ai/i+a 2 /2 J V 



AJ1+AJ2 



( cos(,3) \ ( Ah-AI 2 ^ 2 



(37) 



It can be seen that the support for the distribution W\ 
is the rectangle 

\U\ < C T ° SW - , |V| < ;° S(/3) t , (38) 



where [/ = (A7 X + A7 2 )/e and V = (A7j - A7 2 )/e. 
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FIG. 2. Fidelity decay for the rectangular billiard for vari- 
ous values of the parameter (3=0 (full line), 0.232 (dashed), 
0.3 (dot-dashed) and ty/2 (dotted). The oc 1/t 2 decay is shown 
as a thin dotted line. 

In Fig. 2 we show the decay of fidelity for this system 
for various values of (3. The parameters of the system 
have been chosen as follows: a\ = {-Jh + l)/2, a 2 = 1. 
In all cases the initial phase space density is a hyper- 
rectangle centered around 7l = 1, 7 2 = 1,6 1 = 1 and 
2 = 1 with all sides of length vi x = vi x = vq 1 = v® 2 — 
0.02. The perturbation parameter is e = 3 x 10~ 4 and 
the number of trajectories N = 10 5 . For the above ini- 
tial conditions and parameters ai, a 2 , the critical value 
/3 C which separates the power law and the ballistic fi- 
delity decay is determined by the equality (17) in the 
direction of the U variable. Indeed, when (3 increases, 
— [3fi/dI] -1 (Afi/e) goes outside the support of W\ at 
first along this direction. This gives 



AT 

e 



u 



COs(/3 e ) 

aih + a 2 I 2 



on 
~bT 



An 



(39) 



J u 



where the right-hand side is the U -component of the vec- 
tor 



Therefore we get 

tan(/3 c 







a: 



sm(f3 c )I 2 
sin(/3 c )7i 



aia 2 



(aih + a 2 7 2 ) 



2 ' 



(40) 



(41) 



ai7i + a 2 7 2 



ai7i — a 2 7 2 



Substituting the chosen values of 7i,7 2 ,ai, and a 2 , we 
find that the critical value is equal to (3 C w 0.232. This 
theoretical expectation is confirmed by the numerical 
data of Fig. 2, which show a crossover from a power law fi- 
delity decay (for (3 < (3 C ) to a ballistic decay (for (3 > (3 C ). 
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The results are very similar to the case of the twist map, 
including the fact that, close to the critical value (3 = j3 c , 
the decay is postponed to longer times. Is should be 
stressed that the algebraic fidelity decay, differently from 
the twist map case, is now inversely proportional to the 
square of the time, in agreement with our theoretical ex- 
pectation for a two-dimensional system (see Eq. (14)). 
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FIG. 3. Fidelity decay for the kicked rotor map with 
K = 0.3, e = K' - K = 1(T 5 , and N = 10 4 trajectories, 
using two different initial phase space densities, one centered 
at the point (6 = ir, / = 0.2) (full line) and the other centered 
at (9 = tt,I — 1.2) (dashed line). In the first case the fidelity 
decays ballistically, in the latter case inversely proportional 
to time. Both initial densities are 0.02 wide in the and I 
directions. The oc 1/t curve is shown as a thin dotted line. 

As a last numerical example, we consider the kicked 
rotor map that is given by 

I t+1 =I t + Ksm(e t ), (42) 



IV. CONCLUSIONS AND OUTLOOK 

In this work we have studied the decay of the fidelity of 
classical motion for integrable systems. Our main result 
is the following: for small enough perturbations, the type 
of the decay of fidelity for integrable systems depends not 
on the strength of the perturbation but on its shape and 
on initial conditions. More precisely, the fidelity exhibits 
two completely different behaviors, namely an algebraic 
decay if the perturbation mainly affects the shape of the 
tori, and a faster, ballistic decay, if the main effect of the 
perturbation is to change the frequencies of the tori. Wc 
have also given clear numerical demonstrations of the 
transition between the two types of behaviors, induced 
by changing the shape of the perturbation or the initial 
conditions. 

This result poses interesting questions with respect to 
the quantum mechanical picture. Due to the correspon- 
dence principle, there should exist regimes where both 
types of decay may be observed. It is however expected 
that, for small perturbations, quantum mechanics would 
favor the ballistic type decay, as demonstrated in [11]. 
Indeed the algebraic decay is due to the transitions be- 
tween tori which, for small perturbations, are suppressed 
in quantum mechanics, due to tori quantization and sub- 
sequent gaps between them. The classical-quantum cor- 
respondence will be the topic of further studies. 

This work was supported in part by the EC RTN con- 
tract HPRN-CT-2000-0156, the NSA and ARDA under 
ARO contract No. DAAD19-02-1-0086, the PA INFM 
"Weak chaos: Theory and applications" , and the PRIN 
2002 "Fault tolerance, control and stability in quantum 
information processing" . 



APPENDIX A: ASYMPTOTIC BEHAVIOR 



e t+ i = e t + j t+ i. (43) 

As it is known, for K <C 1 the system is almost inte- 
grable, namely its phase space is dominated by invariant 
tori. There is a stable fixed point at (0 = n, I = 0) and 
a separatrix which divides the phase space into two re- 
gions: a section of librational motion around the stable 
fixed point inside the separatrix, and a section of rota- 
tional motion outside the separatrix. We perturb the 
system by varying K —* K' = K + e. The important 
point is that the type of the perturbation chosen strongly 
affects the frequencies of the tori in the librational sec- 
tion, while it mainly perturbs the shape of the rotational 
tori. Therefore the same system and perturbation should 
lead to two completely different types of fidelity decays, 
power law or ballistic, depending on the choice of the ini- 
tial conditions. Fig. 3 confirms this expectation: if the 
initial density po is inside the separatrix the fidelity decay 
is ballistic, otherwise it is power law. 



The results of the previous Sections do not tackle the 
asymptotic decay of fidelity for integrable systems [23]. 
Indeed, we neglected the contributions to the evolution 
kernel (14) that stem from the fact that the angle vari- 
ables are cyclic. This means that, after a time which is 
oc 27r/e, we need to take into account the contributions to 
(14) not only at ®2t — ©o but also at all © 2 t — @o + 27rk, 
where k is a vector of integer numbers. 

We limit ourselves to the case of a single torus. Of 
course the fidelity f(t) is strictly zero for a single torus, 
and therefore it should be understood that we take the 
limits 

e — > 0, v\ — > 0, with — = constant <C 1. (Al) 
v\ 

Let us consider the initial density po(®) to be defined 
on the whole space (without 2tt periodicity), while 
the kernel Ki(&;t) is defined as the periodic function 
obtained from the original kernel: 



G 



^i(0;i)=^Pi(0-27rk;t) (A2) 

k 

Assuming that the initial density is square normalized, 
the fidelity can be written as 

f(t) = J d n &p* (&)p 2t (&) = J cT*p* (*)p 2t (&), (A3) 
where ~ denotes the Fourier transform: 

p(*) = -jL [ d"0p(0)ex P (-i*0). (A4) 

Since 

P2t(&) = J <T&p {e , )K I (& - &;t), (A5) 
in the Fourier picture this becomes 

fc t (*) = A)(*)£i(*;t)- ( A6 ) 

We may write the original kernel (14) in the simplified 
form 

Pj(®;t) = l/t n Pl (®/t + T), (A7) 

where pi(y) = rWi (- [f^ and T = 

An. The Fourier transform of the kernel (A2) is there- 
fore given by 

ATi(*;t) = £pi(i*)exp(i$(rt - 2nk)). (A8) 

k 

Then the formula 

^ exp(-i27r*k) = ^<J(*-j) (A9) 

k j 

leads to 

#i(*;t) = £>(tj) cxp(ztrj) 5(# - j). (A10) 
j 

This result, coupled with equations (A3) and (A6), finally 
leads to 

/(*) - £ \Po(i)\ 2 M® exp(itrj). (All) 
j 

As we can see, the behavior of fidelity in the limit t — > oo 
is given by the tails of the Fourier transform of the ker- 
nel p\. The origin of the kernel is the projection of the 
perturbed tori onto unperturbed ones, and we expect sin- 
gularities in such a projection. These singularities induce 
a power law decay in the tails of the Fourier transform 
of the kernel and thus are responsible for the asymptotic 
power law decay of fidelity. 

In the single degree of freedom situation, the typi- 
cal singularity of projection to be encountered leads to 
Pi(y) oc \y — yo|~ 1//2 , as it can also be seen in the twist 



map example (22). This type of singularity leads to the 
Fourier transform 

p/($) oc $~ 1/2 exp(-i$y ). (A12) 

Such an expression leads to following asymptotic fidelity 
decay: 

/(*)-/(«>) = 
*- 1/2 £|p°tf)| 2 .r 1/2 exp(*7(ri-2/o)) =t-V 2 z(P), 

(A13) 

where (3 = —Tt + yo and z is some periodic function 
with period 2tt. We note that Eq. (A13) gives an over- 
all cx i -1 / 2 fidelity decay together with a superimposed 
oscillatory behavior. This is the typical asymptotic re- 
laxation of fidelity for a single torus in integrable systems 
with a single degree of freedom. If one considers a finite 
interval of actions Vj, the decay (A13) must be averaged 
over vj, and therefore, due to the oscillatory nature of 
Eq. (A13), it can be faster that i -1 / 2 . The extension to 
the many-dimensional case requires a complex analysis 
of the singularities encountered in the projection of the 
perturbed tori onto the unperturbed ones and is beyond 
the scope of the present paper. 



[1] G. Casati, B.V. Chirikov, I. Guarneri, and D.L. Shep- 
elyansky, Phys. Rev. Lett. 56, 2437 (1986). 

[2] A. Peres, Phys. Rev. A 30, 1610 (1984). 

[3] R.A. Jalabert and H.M. Pastawski, Phys. Rev. Lett, 86, 
2490 (2001). 

[4] F.M. Cucchietti, H.M. Pastawski, and D.A. Wisniacki, 

Phys. Rev. E 65, 045206(R) (2002). 
[5] Ph. Jacquod, P.G. Silvestrov, and C.W.J. Beenakker, 

Phys. Rev. E 64, 055203(R) (2001), 
[6] N.R. Cerruti and S. Tomsovic, Phys. Rev. Lett. 88, 

054103 (2002). 

[7] V.V. Flambaum and F.M. Izrailev, Phys. Rev. E 64, 
036220 (2001). 

[8] Z.P. Karkuszewski, C. Jarzynski, and W.H. Zurek, Phys. 

Rev. Lett. 89, 170405 (2002). 
[9] G. Benenti and G. Casati, Phys. Rev. E 65, 066205 

(2002). 

[10] D.A. Wisniacki and D. Cohen, Phys. Rev. E 66, 046209 
(2002). 

[11] T. Prosen and M. Znidaric, J. Phys. A 35, 1455 (2002). 
[12] Y.S. Weinstein, S. Lloyd, and C. Tsallis, Phys. Rev. Lett. 

89, 214101 (2002). 
[13] W. Wang and B. Li, Phys. Rev. E 66, 056208 (2002). 
[14] Y. Adamov, I.V. Gornyi, and A.D. Mirlin, preprint cond- 

mat/0212065. 



7 



[15] I. Garcia-Mata, M. Saraceno, and M.E. Spina, preprint 

nlin.CD/0301025. 
[16] J. Vanicek and E.J. Heller, preprint quant-ph/0302192. 
[17] T. Prosen and M. Znidaric, J. Phys. A 34, L681 (2001). 
[18] G. Benenti, G. Casati, S. Montangero, and D.L. Shep- 

elyansky, Phys. Rev. Lett. 87, 227901 (2001). 
[19] J. Emerson, Y.S. Weinstein, S. Lloyd, and D.G. Cory, 

Phys. Rev. Lett. 89, 284102 (2002). 
[20] B. Eckhardt, J. Phys. A 36, 371 (2003) 



[21] Ph. Jacquod, I. Adagideli, and C.W.J. Beenaker, Euro- 
phys. Lett. 61, 729 (2003). 

[22] H.C. Corben and P. Stehle, Classical Mechanics (2nd edi- 
tion), (Robert E. Kreiger Publishing, New York, 1974) p. 
244 

[23] For chaotic systems this problem has been addressed 
by G. Benenti, G. Casati, and G. Veble, preprint 
nlin.CD/0208003, to be published in Phys. Rev. E. 



8 



